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1.  Introduction 


Reasonable  understanding  of  the  origin  and  of  the  mathematical,  physical,  and 
practical  roots  of  the  Phenomenological  Mechanochemistry  of  Damage  (PMD)  is 
impossible  without  reading  the  report  on  the  early  history  of  the  PMD.1  Early 
publications  are  deliberately  simplified  to  make  transparent  the  main  object  and 
drama  of  the  PMD:  it  is  the  competition  between  the  (macroscale)  elastic 
accumulated  energy,  from  one  hand,  and  the  (nanoscale)  energy,  accumulated  in 
the  chemical  bonds,  on  the  other  hand.  Because  the  “multiscale  approach”  is  the 
buzz  word  these  days,  I  would  like  to  emphasize  that  the  PMD  theory  is  essentially 
macroscopic  and  phenomenological  approach,  when  treating  both  scales.  The  PMD 
is  the  symbiosis  of  the  damage  theory  as  presented  by  Kachanov2,  and  of  Gibbs’3 
variational  scheme  in  thermodynamics  of  heterogeneous  systems  as  interpreted  in 
the  monograph.4 

Briefly  speaking,  the  PMD  theory  was  bom  out  of  the  practical  necessity  to  suggest 
macroscopic  models  describing  the  widespread  paradoxical  radial  pattern  of 
damage  often  observed  in  terminal  ballistics.  The  radial  cracking  paradox  in 
ballistics  could  be  formally  and  technically  treated  as  the  manifestation  of  the 
typical  Stress  Driven  Rearrangement  Instabilities  (SDRI)  of  phase  interfaces.  The 
SDRI  mechanisms  originated  as  the  natural — both  physically  and 
mathematically — development  of  the  logically  rigorous  variational  scheme  of 
Gibbs.3  It  still  remains  not  understood  by  the  many  followers  of  Gibbs,3  that,  when 
using  the  Gibbs’3  methodology  even  in  the  most  rigorous  and  consistent  way,  we 
arrive  in  99  out  of  100  not  to  novel  physical  phenomena  but  to  novel 
physico/mathematical  paradoxes. 

Discovery  of  paradoxes  is  extremely  important  for  the  progress  of  science  in  many 
respects.  Philosophically  speaking,  we  are  dealing  with  the  apophatic  pattern  of  the 
scientific  study. 

First  and  foremost,  the  SDRI  is  a  set  of  thermodynamic  paradoxes.  Not  more  than 
that — but  not  less  than  that  either.  At  the  same  time,  the  discovery  of  the  SDRI 
paradoxes  demanded  essential  further  developments  of  powerful  technical 
languages  of  compatibility  conditions  of  singular  fields  as  well  as  its 
implementation  into  multidimensional  variational  calculus.  To  the  best  of  this 
author’s  knowledge,  still,  35  years  after  its  appearance,  these  techniques  have  been 
mastered  by  less  than  5  individuals  worldwide.  In  fact,  there  is  nothing 
paradoxical — the  real  pace  of  the  progress  with  fundamentals  is  always  extremely 
slow. 
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Despite  the  fact  that  SDRI  of  phase  interfaces  is  a  set  of  physical  artifacts,  the 
associated  technical  tools  can  and  should  be  used  in  different  applied  areas  of 
physics,  mathematics,  and  mechanics  (but,  of  course,  not  in  the  context  of  the  SDRI 
of  phase  interfaces  in  crystalline  materials).  In  particular,  numerical  modelers  and 
computational  physicists  noticed  in  early  1990s  that  the  SDRI,  in  their  late  stage, 
lead  to  the  appearance  of  cracks.  That  observation  seemed  promising  for 
developing  a  novel  approach  to  formulating  practically  useful  models  of  cracking. 

In  the  terminal  ballistics  and  many  other  potential  applications,  to  which  this  author, 
has  no  direct  relation,  the  mathematics  of  the  SDRI  can  be  immediately  applied  to 
reproduce  the  patterns  typical  for  the  paradox  of  radial  cracking.  However,  one  has 
clearly  distinguish  the  mathematics  of  the  SDRI  (which  is  absolutely  correct!)  from 
the  physics  of  the  SDRI  of  phase  interfaces  (which  is  just  a  physical  artifact!).  Each 
publication,  relating  to  the  SDRI-based  approach  to  radial  cracking,  not 
emphasizing  clearly  and  explicitly  this  remarkable  discrepancy,  would  indoctrinate 
practitioners  with  conceptually  wrong  and  misleading  vision.  It  is  really  unwise  to 
erect  an  important  technical  approach  using  a  theoretical  artifact  as  its  basis. 

I  think,  the  most  productive  strategy  would  be  the  following.  The  SDRI  should 
preserve  its  status  of  a  fundamental  paradox.  Then,  theorists  will  be  prompted  to 
resolve  this  paradox — that  is  the  primary  practical  role  of  any  paradox.  At  the  same 
time,  keeping  the  productivity  of  the  SDRI-based  vision  in  mind,  we  have  to  try  to 
modify  somewhat  both  the  physical  model  and  the  mathematical  technique  of  the 
SDRI.  The  PMD  theory  is  one  of  the  possible  outcomes  of  such  an  approach.  The 
PMD  does  not  mention  explicitly  or  implicitly  any  phase-transformations 
whatsoever.  It  does  not  mention  explicitly  the  SDRI  whatsoever.  It  does  not  use 
directly  any  techniques  associated  with  the  SDRI  theory  and  even  with  the  Gibbs’3 
paradigm.  Yet,  it  uses  the  clear  variational  paradigm,  typical  for  the  classical  statics. 
Instead  of  the  SDRI,  it  uses  a  more  traditional  concept  of  thermodynamic 
instability.  It  allows  to  reproduce  the  radial  cracking  pattern. 

Keeping  this  promising  progress  in  mind,  we  might  begin  thinking  of  enriching  the 
original  simple  PMD  model  with  additional  physical  features.  In  this  report,  we 
formulate  the  extension  of  the  PMD  theory  allowing  one  to  include  ponderomotive 
electrostatic  or  magnetostatic  forces  into  the  general  scheme. 

In  this  report,  we  formulate  the  extension  of  the  PMD  theory  allowing  one  to 
include  ponderomotive  electromagnetic  forces  in  the  general  scheme.  Our  approach 
of  including  the  ponderomotive  forces  is  based  on  earlier  publications.5  7  One  of 
our  ultimate  goals,  as  we  foresee  it  now,  is  to  suggest  a  novel  model  of  electric 
breakdown8,9  based  on  the  extended  PMD  approach. 
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There  are  thousands  of  publications  devoted  to  analysis  of  ponderomotive  forces 
including  classical  textbooks10-13  among  others.  The  specific  feature  of  those 
publications  is  the  analysis  of  the  problem  in  the  framework  of  linear  elasticity. 
Toupin14  was  one  who  started  the  analysis  in  the  full  nonlinear  framework.  In  our 
opinion,  he  made  a  strategical  mistake  though,  relying  on  usage  of  the  Lagrangian 
framework.  Thus,  he  considerably  narrowed  the  applicability  of  his  theory  since  the 
domain,  occupied  by  the  substance,  typically  has  a  negligible  size  as  compared  with 
the  whole  space  occupied  by  the  electromagnetic  field.  Unfortunately,  his  mistake 
was  made  by  many  of  his  followers. 

The  approach  that  is  presented  in  the  work  found  in  references  5-7  and  15  relies  on 
the  usage  of  the  Eulerian  description  and  it  is  free  of  the  drawbacks  mentioned  in 
the  Lagrangian  approach  of  Toupin  and  his  followers. 

2.  Quasi-Static  Master  System  for  Deformable  Polarizable 
Solids 


Consider  a  deformable  polarizable  substance  in  a  thermostat  maintained  at  fixed 
temperature  T° .  In  what  follows  the  parameter  T°  will  be  omitted  from  all  the 
relationships.  The  free  energy  y/  per  unit  mass  is  given  by  the  following  formula 

y/  =  y/(ViUj,Pk,KM),  (1) 

where  Pk — the  polarization  vector,  Ui — the  displacement  vector,  km — the 
damage  parameters. 

This  substance  reacts  on  the  external  load  by  generating  elastic  strain  and 

electric  polarization  Pk .  Under  the  action  of  these  agents  the  integrity  bonds  of  the 
substance  gets  broken.  Since  there  are  different  bonds  in  solids  it  might  be 
worthwhile  to  use  different  damage  parameters.  In  the  following  though,  we  will 
use  for  brevity  only  1  damage  parameter  k  . 

We  assume  that  the  mechanical  and  electrostatic  equilibrium  in  the  system 
establishes  much  faster  as  compared  with  the  time-scale  of  establishing  “chemical” 
equilibrium.  As  before,  by  chemical  equilibrium  we  understand  the  equilibrium 
with  respect  to  variation  in  the  magnitude  of  the  damage  parameter  k  . 

We  suggest  analyzing  a  slow  evolution  of  this  system  with  the  help  of  the  following 
master  system.  This  system  includes  the  system  of  electrostatics,  consisting  of  the 
bulk  equations: 

Vp  =  0  (2) 
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and 


zWcVjEk=  0,  (3) 

where  E'  and  D'  are  the  electrostatic  field  and  displacement,  satisfying  the  identity 

D‘  -  E‘  +  4nPl .  (4) 

At  the  internal  and  external  boundaries  these  vectors  should  satisfy  the  following 
boundary  conditions: 

|>'lfl=0  (5) 


and 

[D']>,=0,  (6) 

where  Qi  and  A,  are  the  tangent  and  nonnal  vectors  to  the  boundary. 


The  bulk  equations  of  the  electrostatic  and  mechanical  equilibrium  read 


dys  77 
P — -  =  E. 

dP‘ 


(7) 


and 


V.  Kmk  =  0 


(8) 


where  the  Aleph  tensor  is  defined  as  follows: 


=  p 


dy/ 


av  u. 

m  j 


(spvp) 


—  z 


mk 


1 


1 


— E.D1 - E.E' 

\4  n  8  n  j 


\  1 

+ — DmEk 
4  n 


(9) 


(see  references  5-7  and  15). 

By  definition,  at  the  coherent  interface  the  displacement  vector  should  be 
continuous 

(10) 


together  with  normal  components  of  the  Aleph  tensor 

[«”*]>„  =0.  (11) 


The  same  condition  (Eq.  11)  should  be  satisfied  at  the  interface  substance/vacuum 
(obviously,  the  free  energy  of  the  vacuum  should  be  ignored). 
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At  last,  we  postulate  the  following  equation  for  the  rate  of  damage: 


dt 


dt  1  die 


(12) 


where  V  (z,t)  is  the  velocity  field  of  substance  satisfying  the  identity 


(13) 


Here  K  is  the  rate  of  kinetics  state  function.  It  should  be  chosen  based  on  existing 
experimental  data.  The  only  universal  constraint  on  it  is  the  positivity:  K  >  0 .  Also, 
from  comparison  with  experimental  data  one  has  to  choose  the  free-energy  density 
function  i//  =^(V  JJ pPk ,ie).  Of  course,  it  is  natural  to  demand  that  the  system 

Equations  1-13  obeys  the  laws  of  thermodynamics.  This  demands  imposing  some 
constraints  on  the  admissible  functions  y/  =  ^(V  JJ .,Pk ,k).  These  constraints  are 

known  as  thermodynamic  inequalities.  Still,  a  lot  of  freedom  remains  in  the  choice 
of  y/(y  JJ j,Pk ,k).  Thermodynamics  is  helpless  in  generating  further  constraints. 

3.  The  PMD  of  Nonpolarizable  Deformable  Solids 

For  this  special  case,  the  energy  density  does  not  depend  upon  the  polarization 
vector 


(//  =  ¥(VUpk) 


(14) 


the  general  master  system  reduces  to  the  following: 

=  0, 


(15) 


where  the  Aleph  tensor  is",k  reduces  to  the  following  form: 


(16) 


which  is  nothing  but  the  Cauchy  stress  tensor  in  the  Eulerian  variables  (see 
reference  4  and  the  references  therein). 

By  definition,  at  the  coherent  interface  the  displacement  vector  should  be 
continuous 


(17) 


5 


together  with  normal  components  of  the  Aleph  tensor 


[tr‘]X=o. 


(18) 


The  same  condition  (Eq.  11)  should  be  satisfied  at  the  interface  substance  vacuum 
(the  free  energy  of  the  vacuum  should  be  ignored). 


And  finally,  we  postulate  the  following  equation  for  the  rate  of  damage: 

d/c(z,t)  dy/(VU  nK) 

— \J-L  +  v,WiK  =  -K  J  , 

dt  die 

where  V  (z,t)  is  the  velocity  field  of  substance  satisfying  the  identity 


(19) 


dU’ (z,t)  ./  .  .  \ 

- ^-L  =  V](d\-V ,Ul)  . 

dt  \  j  j  J 


(20) 


Thus,  in  addition  to  the  exact  nonlinear  master  system  of  the  PMD  in  the 
Lagrangian  framework16-18,  we  attained  the  exact  nonlinear  master  system  of  the 
PMD  in  the  Eulerian  framework. 


4.  Quasi-Static  Master  System  for  Rigid  Polarizable  Solids 

The  above  system  remains  meaningful  even  in  the  case  of  rigid  (i.e., 
nondeformable)  substances.  In  this  case,  the  general  master  system  for  deformable 
polarizable  substances  reduces  to  the  following: 

,//=,//(P\0  (21) 

V,.D'  =  0  (22) 

zjkVjEk=  0,  (23) 

where  E'  and  D‘  are  the  electrostatic  field  and  displacement,  satisfying  the  identity 

D‘  -  E‘  +  4ttP‘  .  (24) 

At  the  internal  and  external  boundaries  these  vectors  should  satisfy  the  following 
boundary  conditions: 

[>'!&=  0  (25) 

[D']"A,.=0,  (26) 

where  Qf  and  Nt  are  the  tangent  and  normal  vectors  to  the  boundary. 
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The  bulk  equations  of  the  electrostatic  and  mechanical  equilibrium  read 

r 

p-L-  =  E.. 

8F 

Also,  we  postulate  the  following  equation  for  the  rate  of  damage: 

dK(z,t)  _  K8y/{Pk,K ) 

8t  d/c 

In  this  case,  the  density  p  in  Eq.  27  is  just  an  unchangeable  constant. 


(27) 


(28) 


5.  The  PMD  for  the  Linear  Isotropic  Substances 


Consider  the  simplest  free  energy  density  of  rigid  dielectrics 

AF‘^)=^yrpEp‘*c(K). 


(29) 


Then,  thermodynamics  imply 


D‘  (Pk ,tc)  =  E‘  +  AkP'  =  A7l£[K^  P<  =£(K)Ei 
V  '  £(/r)-l 


(30) 


The  electrochemical  master  system  reads 


8t 


±-E’E,e.(K)+C,(K)\  , 


(31) 


and 


V,(££')  =  0  . 


(32) 


Let  us  choose  the  simplest  model 


s[k)  =  \  +  s°  (\-ax),  C(/r)  =  y(/r-K*°)  . 


(33) 


Then,  the  system  can  be  rewritten  as 


1  d/c  „  1  / 

- =  sa — E'E  .-y\k-k 

K  8t  8  n  1  K  ’ 


(34) 
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and 


V;.{(l  +  £°(l-aK-))£'}  =  0. 


(35) 


6.  Spatially  Uniform  Equilibrium  Configuration 


implying 

k*  =  k°  ^  —  E^E' 

4  nx 

and 

s*  =  \  +  s°  —E*JE* . 

\  !  4^, 


(36) 


(37) 


(38) 


For  negligibly  small  initial  damage  k  we  get,  respectively,  the  following  values  of 
the  damage  and  electric  permittivity 


and 


* 


K 


£  a 


E*JE* 


o2  2 

„  £  a  , 

=  l  +  £ - E  'E 


4  nx 


J 


(39) 


(40) 


7.  System  of  Small  Disturbances  in  Vicinity  of  Spatially 
Uniform  Configuration 

We  obtain  the  system  of  small  disturbances  in  vicinity  of  the  unifonn  equilibrium 
configuration  by  linearizing  the  bulk  master  system  (Eqs.  3 1  and  32).  We  get 

— —  =  £°a  —  E*'E.  -  xi< ,  (41) 

K  dt  An  1 

and 

Vi(-aice°E*i  +  £*&)  =  0.  (42) 
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In  1 -dimensional  case  the  system  (Eqs.  41  and  42)  reads,  respectively,  as 


and 


1  die 
K  dt 


a  a 
4  K 


E*E  -  xk  . 


(43) 


-aks°E*  +s*E  -  A(t)  . 


(44) 


We  are  looking  for  solutions  of  the  linear  system  (Eqs.  43  and  44)  in  the  form 

ie(x,t)  =  Pe,1,+ihl,  E(x,t)  =  Qe'lt+ila ,  (45) 

where  k  is  the  spatial  wave-vector,  ^  is  the  rate  of  growth  of  small  disturbances, 
whereas  P  and  Q  are  constants. 

Substituting  Eq.  45  in  Eqs.  43  and  44,  we  get 


K 


+  2T 


p^L?Le*q  =  o. 

An 


(46) 


and 


as°E*P-s*Q  =  0.  (47) 

The  linear  system  (Eqs.  47  and  48)  with  respect  to  the  constants  P  and  Q  has 
nonvanishing  solutions  only  when  its  determinant  vanishes.  Equating  the 
detenninant  to  zero  we  arrive  at  the  secular  equation: 

p2  2 

—Et2=  0.  (48) 

An 

In  view  of  Eq.  47,  we  arrive  at  the  following  formula  for  the  rate  ij  of  the  bulk 
disturbances: 


n_ 

K 


-2T  + 


o2  2 

£  CC 

4  ns* 


(49) 


When  the  rate  of  growth  tj  is  less  than  zero  the  disturbances  decay  exponentially. 
When  ij  is  greater  than  zero  they  grow  exponentially  in  time.  Because  of  this 
exponential  growth  the  amplitudes  of  the  disturbances  become  so  big  that  the 
linearization  becomes  invalid.  However,  the  linear  analysis,  on  its  own,  allows  to 
conclude  that  the  uniform  configuration  becomes  unstable. 
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The  formula  in  Eq.  49  clearly  demonstrates  the  competition  in  the  influences 
between  the  chemistry  of  the  bond’s  damage  and  ponderomotive  forces.  The 
chemistry  of  the  bond’s  damage  produces  the  stabilizing  effects,  whereas  the 
ponderomotive  forces  try  to  destabilize  the  system  by  enhancing  the  damage.  We 
call  as  critical  the  electric  field  at  which  the  rate  //  crosses  the  zero  (neutral)  value. 
Using  Eq.  50  we  arrive  at  the  following  formula  of  the  Ecrit 

Elu=X^r{}  +  e°(\-aK°)Y  (50) 

s  a  v  v  ” 


In  the  case  of  the  small  initial  damage  k°  the  last  fonnula  reads 

2k  l  +  £° 


Elit  =  Z: 


a~ 


(51) 


8.  Conclusions 


We  discussed  a  phenomenological  variational  approach  to  the  problems  of 
electrostatics  and  magnetostatics.  Although  we  explicitly  focused  on  electrostatic 
heterogeneous  systems,  almost  all  of  the  presented  results  are  valid  for 
magnetostatics.  Our  approach  is  close  to  the  Gibbs’3  minimum  energy  approach,  as 
it  was  interpreted  in  Grinfeld  1994. 4 

Although  the  suggested  approach  is  based  on  logically  consistent  derivation  from 
the  Gibbs’3  energy  principle,  its  implications  should  not  be  treated  as  reliable 
physical  facts.  Potentially,  the  discrepancy  of  many  of  those  implications  against 
experiments  and  observations  should  be  treated  as  paradoxes.  We  believe  that  it  is 
important  to  remember  that  paradoxes  always  played  an  extremely  important  role 
in  the  development  of  fundamentals  of  physical  sciences. 
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